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Finding Exact Values For Infinite Sums 
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Introduction In the 1995 February issue of Math Horizons I. Fisher posed the following problem 

[1]: 

From the well-known results, 

E — = — and > 

it follows that 
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^ n 2 + n/2 6 
Find £/ie exact value of the convergent sum. 

This paper offers a solution method that allows one to find exact values for a large class of 
convergent series of rational terms. 

In the next section, we first illustrate the method for a special case. We then describe the general 
result pointing out further generalizations of the method and we finally end with a brief discussion. 

A Special Case Consider the series 



S(a,b) = J2 



^ (n + a)(n + b) 

71=1 v ' v ' 

where a ^ b and neither a nor b is a negative integer. Sums of this form arise often in problems 
dealing with Quantum Field Theory (p. 89ff., Ref. [2]). 
Decomposing each term of (|l|) in partial fractions gives 

5 M) = E ( —u - —) ■ (2) 

v ' ' a-b \n + b n + a > 

71=1 

Now we use the identity 



1 

A 

Therefore for a,b > —1 we have 

S(a, b) = - 
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e~ Ax dx , A > . (3) 
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b J a 



dx . (4) 



In deriving the last result, we made use of the monotone convergence theorem (p. 318, Ref. [3]). In 
particular, the integrand in the second line of equation (f|) consists of a non-decreasing sequence of 
non-negative functions and therefore we can swap the indicated operations of taking the limit and 
performing the integration. 

Making the change of variable t = e~ x in the integral of equation (|]) gives a more symmetric 
result: 

1 f 1 t b - t° 



S(a,b) = / -r—dt. (5) 

a-b J 1 - t 



Some comments are in order here: 
• Although the integral 



t y , 
dt 



o 



l-t 

diverges, the integral of equation ([5]) converges for a ^ b and a,b > —1. 



The problem of Fisher that appeared in Math Horizons corresponds to a = 1/2, 6 = 0. In this 
case, we have 

v 1 2 r 

^ n(n+l/2) J i + VT 

After the change of variables t = u 2 , it is easy to calculate the integral: 

— — l —— - = 4 / ( 1 —\ du = 4(1 - In 2) ~ 1.227 . 

When the two numbers a and b differ by an integer k, i.e. a = b + k, then the sum ([!]) 
"telescopes" and it can be easily calculated from (H): 

This can be used as a consistency check of formula (|5|). Indeed 

-1 +k 1 i /-l fc - 1 i k - 1 j-i+a+l 



S(a,a-k) = ijf ^_ifcft=ijf E^^=^E 



z + a + 1 
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in agreement with the last result. 

Now we note that one can express the result (|||) in another equivalent form — namely, using the 
well known representation (see p. 258, Ref. [4]) of the digamma function ip{z): 

1,{z) = ^lnr(z) = - 7 - J^-^^dt, (6) 



where 7 is the Euler-Mascheroni constant and T(z) is the gamma function. This is motivated by the 
fact that part of the integrand in the r.h.s. of equation (||) is similar to the integrand in (^j). In fact, 
we have 

b — a 

There are many useful identities involving the digamma function (p. 258, Ref. [4]). For example 

il>(l + z)=ip(z) + - . (8) 
Moreover, the exact value of ip{z) is known for several values of z: 

V(l) = -7, V(1/2) = -7-21b2. (9) 

Equations (||) and @ can be used to evaluate S(a, b) exactly for many values of a and b. For 
example, if a = 1/2, b = we find 



5(0,1/2) 



2[V(3/2)-V(l)] = 2[^(l/2) + 2-V(l)] 
2[-7-21n2 + 2 + 7 ] = 4(l-ln2), 



in agreement with our previous result. 

When a = b in ([TJ), the sum can still be calculated. We consider two approaches. The first 
approach is to repeat the calculations presented above but observing that the basic equation (H) has 
now to be modified in the form 
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A* 

Following the same reasoning, we find 
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-Ax 



dx 



A > 0. 
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(10) 



Alternatively, we can obtain the same result by taking the limit b — > a in (0): 

1_ f 1 t b -t a -1 



lim 

b— »a 
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dt . 
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Without loss of generality, we can assume that — 1 < a < b. We notice that for < t < 1/2 

t b -t a 

< 2t° 



'IV, 
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while for 1/2 < t < 1 the integrand of (|5|) is bounded; let M be its supremum in this subdomain. 
The function 



g(t) 



g , if < t < 1/2 

m", if 1/2 < t < 1 
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is integrable in [0,1] and therefore the dominated convergence theorem (p. 167, 321, Ref. [3]) can 
be used to interchange the operations of the integral and the limit in ([TT]) : 
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Also, from equation (|7|) we find 



, t b -t a -1 , f 1 t a lnt , 
lim at — — I at . 



.n + a) 2 J n b-*a b — al — t / n 1 — t 

n=l ' 



El dip(z) 
In 4. 



z=a+l 
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The functions 

^\z) = ^rlnrW = -jf dt * /_7 , n = l,2,... (12) 

are known as polygamma functions (p. 260, Ref. [4]). Several identities for the polygamma functions 
are known (p. 258ff., Ref. [4]). For example 

^W(l) = (-l) n+1 n \((n+ 1) , 
^ (n) (l/2) = {-l) n+1 n\ (2 n+1 - l)C(n + 1) , 



where ((z) is the zeta function: 



z n+l 



oo ^ 

C(z) = £ - , Re, > 1 . 

n=l 



As an application of fllPp, we obtain the well known result 



1 hit , 7T 2 



dt = — ~ 1.645 . 
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Formula (|T^) also implies the less well known result 

> --r = - / dt = 3( 2 -4 = 4 = 0.935. 

Z^(n + l/2) 2 J l-t ^> 2 

The General Case After our preceding discussion, we can now establish a more general result. 

Let 

Qjv_ 2 (w) 

where QN-2{ n )-, Pn^) are two polynomials in n of degree N — 2 and A" respectively. We shall assume 
that P/v(ra) is expressible in the form 



Pjv(n) = in + ax) mi (^ + « 2 ) m2 • • • (n + a k ) m " , 
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with all dj, i = 1, 2, . . . , k distinct real numbers none of which is a negative integer. This ensures the 
convergence of S. Then for any polynomial 

QN-2{n) = Cat_ 2 n N ~ 2 + Cat_ 3 n N ~ 3 + . . . + c , 

the sum S is written in terms of partial fractions: 

x k mi 

S(ai, . . . a k ; c , . . . , c N _ 2 ) = 2, / v / > 7 — I ^ m > 

?1=1 1 = 1 J = l 

where the constants Ajj are uniquely determined by the partial fraction decomposition of each sum- 
mand. In particular, notice that since there is no term of degree N — 1 in QN-2{n), 



Ai = . (13) 



1=1 



Using the identity 

1 1 /*°° 

x L - x e~ Ax dx, (14) 



A L (L-1)\J 
we write the series in an integral form valid only if > — 1, Vz: 

OO k TTii . pOG 

a fe ; Co,..., CJV _ 2 ) = EEEt-TTv/ ^ dx ■ ( 15 ) 

n=l j=l j=l ^ 

Working in a similar fashion as in the derivation of equation (HI), we find 



k rrii 



S(a h . . . a k ; c , . . . , c N - 2 ) = (j -^l)! / 



i=l j=2 
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1 - e" x 
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where we have taken extra care for the j — 1 term (by using the condition (|T3"D) in order to guarantee 
the convergence of the corresponding integral. This is our result in an integral form. We can also 
express it in terms of the polygamma functions fll2[): 



k m, , ^ 



S(cn, ...a k ;co, ...,c N - 2 ) = ^ r - IV Aij ^ 1) ( a i + 1 ) ? ( 16 ) 

i=l j=i v /• 

where we have defined tp(°\z) = ip(z). 

As a straightforward application of our method, let us consider the following examples: 

Example 1 : 

+oo 



S( a , -a) = ~ri 



— a 2 



where a is not a positive integer and satisfies the inequality a > — 1. Using the formula ( |16| ) 
we find 

, , ip(a + 1) -i/j(-a + 1) 
S{a,-a) = , a^O. 

This result can be further simplified if we make use of the functional relation 

Z + 1) = 1p(z) + 71 COt(7T2:) , 

in conjunction with (§). Then 

"1 



S(a,-a) = - 



Example 2 : 



— 7i cotlna) 
a 



N 

where N > 2 and a > — 1 



1 (— 1 V v 



Example 3 : 



«(o, 0,3/2) = E^TTT72) 



Using formula (pi]), we find 



TV 2 



Example 4 : 



5(0,0,3/2) = 4^(1) + 2^ (1) (1) -4^(3/2) = 8(l-ln2) = 0.835. 

3 



*(M.l/2)= g (n+1) , (n+1/2) 



Using the formula ([16]) we find 

7T 2 

5(1,1,1/2) = [4^(2) -2^ (1) (2) -4^(3/2)] = 2(41n2-l)-— = 0.255. 

Finally, the reader is invited to write down the values for other infinite sums of the form presented 
above. 

Discussion Before we finish, we would like to point out that the identities ([|) and (|l"4D we used 
in our derivations express the quantities 1/A and 1/A L as the Laplace transforms of 1 and x L ~ l 
respectively. In general, if f(s) is the Laplace transform of g(x), 



f(s) = I e sx g(x)dx 

then the sum 



Si = J2 /(») 



ngT" 
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where JcZ, can be written in the form 



Si = jH g{x) (j>— j dx, 



assuming that the operations of summation and integration are interchangable. Assuming moreover 
that the sum inside the parenthesis can be performed explicitly, we have thus obtained an integral 
representation of Sj. 

The Laplace transform has been proved a very valuable tool in the solution of differential equa- 
tions. Unfortunately, in the summation of series, the Laplace transform does not enjoy the same 
popularity. In this paper, we have tried to present some of the limitless possibilities that the method 
offers. We propose our reader to solve the following problem: 

Problem 

(i) Show that 

(~l) rt+1 = f 1 t a dt 
n + a J 1 + t 

n=l 

(ii) Using the previous result, derive the well known result 

00 ( — l)n+l 

= ln2 = 0.693 , 



n 

n=l 



and the less known result 



E izi W = 2-£ = 0.429. 
^ n+ 1/2 2 

n=l ' 



We hope that this will motivate him/her to explore more aspects of the method presented in this 
paper and establish many additional results. 
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